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While entanglement is believed to be an important ingredient in understanding quan¬ 
tum many-body physics, the complexity of its characterization scales very unfavor¬ 
ably with the size of the system. Finding super-sets of the set of separable states 
that admit a simpler description has proven to be a fruitful approach in the bipartite 
setting. In this paper we discuss a systematic way of characterizing multiparticle en¬ 
tanglement via various relaxations. We furthermore describe an operational witness 
construction arising from such relaxations that is capable of detecting every entangled 
state. Finally, we also derive an analytic upper-bound on the volume of biseparable 
states and show that the volume of the states with a positive partial transpose for any 
split rapidly outgrows this volume. This proves that simple semi-definite relaxations 
in the multiparticle case cannot be an equally good approximation for any scenario. 
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I. INTRODUCTION 


Without a doubt entanglement can be considered one of the most important concepts 
in quantum physics, clearly distinguishing quantum systems from classical ones. It can be 
harnessed to enable novel ways of processing quantum information in numerous ways, from 
communication to computation. Many of these operational tasks require an operational 
detection or even quantification of this indispensable resource. While in bipartite systems of 
low dimensions this can be achieved in an efficient way, the complexity of the characterization 
of entangled states makes a complete and computable framework of entanglement detection 
impossible in high dimensions and thus also for multipartite systems^ 2 . 

A possible way for deriving statements on the presence of entanglement is to discard 
the actual complex structure of the border between separable and entangled states and 
try to find good approximations that admit a more amenable description. In the bipartite 
case positive maps play a central role in such approximations: All separable states remain 
positive semi-definite under application of a positive, yet not completely positive map to one 
of its subsystems. The most well-known example of such a map is the partial transposition: 
This map generally changes the eigenvalues of a matrix, but separable states have a positive 
partial transpose (PPT). The approach of positive maps allows for a characterization of 
super-sets of the set of separable states using techniques from semi-definite programming, 
and it nevertheless captures the whole structure: A state remains positive under all positive 
maps if and only if the state is indeed separable^-. 

In order to gauge the efficiency of certain maps for characterizing entanglement, one of 
the most relevant issues is how the volume of states that remain positive under the map 
in question compares to the volume of separable states^. The sobering and non-surprising 
answer from bipartite systems can be gained from convex geometry considerations and shows 
that for all known maps the states that remain positive are most likely entangled in high 
dimensions^. For small dimensions, however, a given positive map can detect a large fraction 
of entangled states. This is, for instance, true for the partial transposition, which delivers a 
necessary and sufficient criterion for 2 x 2 and 2x3 systems. 

In multipartite systems the characterization of entanglement constitutes an even greater 
challenge. Since partial separability of multipartite states can no longer be defined as a 
purely bipartite concept, the application of positive maps to subsystems alone can reveal 
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little more than entanglement across a fixed partition of the multipartite state. Nevertheless, 
recently several works succeeded in defining suitable mixtures of positive maps, which can 
be used to develop strong criteria for genuine multiparticle entanglement^. 

In this paper we first develop a framework that allows for the semi-definite relaxation of 
partially separable states, opening the possibility for harnessing well developed techniques 
based on positive maps to detect genuine multipartite entanglement (an approach which 
has already been shown to yield useful results with different relaxations in Ref.). We 
achieve this goal by first formally defining semi-definite relaxations of partially separable 
states using positive maps. 

Due to the formulation as semi-definite programs these constructions yield versatile cri¬ 
teria for detecting multipartite entanglement in low dimensional systems. To unlock these 
powerful techniques for more complex quantum states we proceed to discuss a recently intro¬ 
duced program of lifting bipartite witnesses^. We prove that it is always possible to exploit 
witnesses that only reveal bipartite entanglement in order to construct witnesses for genuine 
multipartite entanglement. This facilitates this notoriously hard problem, and we showcase 
this technique with some exemplary multipartite entangled states. 

In a second step, we ask which fraction of genuinely multipartite entangled states can 
be detected with such relaxation methods. We prove an upper-bound on the volume of the 
biseparable states, and a lower-bound on the volume of a set of states that can never be 
detected with relaxation methods based on the partial transposition. For large dimensions, 
both values deviate significantly. This shows that while the relaxation approaches are strong 
for small systems, they fail to deliver a good approximation in the general case. 


II. CHARACTERIZING RELAXATIONS OF SEPARABILITY WITH 
SEMIDEFINITE PROGRAMS 

The most straightforward relaxation of separability in multipartite systems is again given 
by positive maps. Trying to justify this assertion is the object of the current section. 

Let us start with some basic definitions and notations (see e.g. Ref.-— for a general review 
on these notions). On a multipartite system, given a bipartition {b\b} of the subsystems, we 
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denote by Sb the set of states which are biseparable across this cut, i.e. 

P e Sb if p = X] ®106°) <06°| ® l4 0 H4 0 l> (!) 

i 

where {(&}* is a convex combination, and for each i, |0^), |</>^) are pure states on the 
subsystems in b , b respectively. We then define the set «S( 2 ) of biseparable states as being the 
convex hull of {<S 6 } b , i.e. 

p G 5 ( 2) if P = (2) 

b 

where {pb}b is a convex combination, and for each 6, cr fe e If a state is not biseparable it 
is called genuinely multipartite entangled GME. 

This definition admits a simple relaxation with a positive-semidehnite characterization. 
Given, for each bipartition {b\b}, a positive map A&, acting on the substystems in b, we 
define the set Tl{\ b } b of {Ab}b -relaxation of <S( 2 ) by 

P e ^{a 6 } 6 if P = ^2PbVA b , (3) 

b 

where {pb}b is a convex combination and for each 6, At <8) 1 &[ctaJ ^ 0. 

Such a relaxation carries the operational advantage that it can be approached via semi- 
definite programming (SDP). Besides, this definition can easily be extended to fc-separable 
states by applying different maps to the induced partitions. As we are however mainly 
interested in characterizing the strongest form of multipartite entanglement we focus here 
on the distinction between biseparable states and genuinely multipartite entangled ones. 

For instance, these relaxations can be particularly useful when optimizing convex func¬ 
tions over the set of biseparable states. Indeed, for any function /, we trivially have that, 
for any set of positive maps {A b}b, 

min /(er) ^ min /(cr). (4) 

o-€S (2) a&TZ {Ab}b 

In particular, for any function / satisfying f(\pi + (1 — A)p 2 ) ^ A/(pi) + (1 — A)/(p 2 ), 
Eq. © provides a relaxation of the optimization of / over Sm which can be cast as an SDP. 
One of the most straightforward applications of such a strategy is to testing whether a given 
density matrix p is indeed biseparable, i.e. applying it e.g. to / : a i —> (Tr [(p — <r) 2 ]) 1 ' / “. It 
yields in that case the equivalence 

V {A b } b , min Tr [a(a - 2p)l ^ Tr [p 2 ] «=> p6 5 (2 ). (5) 
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While this defines, in theory, a necessary and sufficient program for deciding whether a 
given state is biseparable, checking all possible sets of positive maps is of course not feasible. 
However even making one particular choice, such as the transposition map for instance, has 
already proven to yield very strong witnesses, and through the dual of the program one 
can additionally often extract analytical constructions for important classes of states^. To 
test the further prospects of such SDP approaches we have first programmed a Choi-map- 12 
relaxation test in the following way: 

max 6 ’ s.t. (Xi ^ si, <T 2 ^ si, CX 3 = P — (J\ — <J 2 , of 1 ^ si, erf 2 ^ si, a 3 3 ^ ) ( 6 ) 

where C denotes the Choi map. While not being quantitative as the above program it can be 
easily implemented in MATLAB, using the packages YALMIP— and the SEDUMI solver—. 
If an s ^ 0 is found, the program also returns the decomposition elements a* and can thus be 
used to End actual Choi-positive decompositions. As a first test of such a program we have 
chosen the family of states p({A a }) introduced in Ref.-. Setting all X a = X it was shown 
in Ref.- that these states are GME for 0 < A < f. Testing the program on this family of 
states returns a Choi-positive decomposition for | A < 1, thus showing that on the set 
of Choi-positive mixtures the witness presented in Ref.- is weakly optimal. It furthermore 
proves that for A ^ A < 1 the state can be decomposed into states which are either PPT or 
Choi-positive. 

We have also extended the program to demand simultaneous positivity under partial trans¬ 
position and the Choi-map. That is formally, we looked for 

maxs s.t. (T 1 ^ si, cr 2 ^ si, cr 3 = p — <7i — cr 2 , of 1 ^ si, a f 2 ^ si, erf 3 ^ si, 

a\ x ^ si, erj 2 ^ si, a^ 3 ^ si. (7) 

With this program it is easy to check that indeed every state between 0 < A < 1 is genuinely 
multipartite entangled, which proves the versatility of this approach. 

III. CONSTRUCTING MULTIPARTICLE WITNESSES FROM 
BIPARTITE WITNESSES 

While the semidehnite program presented before technically gives sufficient criteria for 
deciding biseparability, it becomes quickly intractable beyond a few qubits. There are, 
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however, frequent situations in which one can use some additional knowledge to facilitate 
witness constructions. Genuinely multipartite entangled states of course also have to be 
entangled across every bipartition of the system. And since the construction of bipartite 
entanglement witnesses can be a rather straightforward affair (e.g. through positive maps) 
one can ask whether there is a possibility to construct multipartite entanglement witnesses 
from a collection of bipartite operators. 

A. A systematic construction for lifting bipartite witnesses 

In Ref.- the authors introduced a general witness construction method, which enables the 
construction of multipartite entanglement witnesses from a set of bipartite witnesses across 
every possible bipartition. Such a problem can be formalized as follows. 

Given, for each bipartition {b\b}, a witness W b: i.e. a self-adjoint operator with the 
property that Tr (W b a b ) ^ 0 for all a b G S b , we are looking for an operator Wgme with the 
following property 


V 6, W G me > W b , (8) 

where, for any self-adjoint operators A , B, we say that A ^ B \i A — 5 is a positive 
semidefinite operator. 

The construction of Ref.- constitutes one analytic answer to this question in the following 
way. Assume a common operator Q for every bipartite witness W b , i.e. 

W b = Q + T b . (9) 

Note that the only obstruction to the existence of such a non-trivial operator Q is that the 
operators {W b } b have the property that for each b there exists b' such that the non-zero 
eigenspace of W b is contained in the zero eigenspace of Wy. Otherwise, one can write down 
a formal solution of the above problem as 

Wgme= Q + £pi] + , (10) 

b 

where, for any self-adjoint operator A, we denote by [A] + the projection of A onto the 
positive semidefinite cone. It is now easy to see that condition (jHJ) holds for every bipartite 
witness W b . 
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Of course, the crucial issues here are first the one of the optimal choice of Q and second 
the one of the generality of such a result. Using a special choice for Q, the authors of Ref.— 
show that there exist genuinely multipartite entangled states and a set of bipartite witnesses 
for which the expectation value of Wqme is negative, proving that this construction indeed 
succeeds in generating witnesses for multipartite entanglement. Here we prove first of all 
that, for every genuinely multipartite entangled state Pgme , there exist a set of optimal 
bipartite witnesses for which this construction yields a multipartite entanglement witness 
Wgme such that Tr (pgmeWgme) < 0. 


Theorem III. 1 For every genuinely multipartite entangled state Pgme, there exists a set 
of weakly optimal bipartite entanglement witnesses {Wb}b such that 


Tr ^ pgme f <3 + l Tb \ +11 


< 0 , 


( 11 ) 


where the operators Q and {Tb}b are defined by equation 


Proof: Let Q be a weakly optimal witness for multipartite entanglement, which means that 
V cr e <S( 2 ), Tr (Qcr) ^ 0, and 3 b 0 , 3 p bo E Sb 0 : Tr(Qpb 0 ) = 0. It is sufficient to show that 
such a witness can always be constructed from bipartite weakly optimal witnesses, as then 
clearly every multipartite entangled state can be detected by the construction. 

From the assumption it follows that 


V b, min Tr (Qpb) ='■ otb IP 0. 

Pb£$b 


( 12 ) 


As this minimization is convex in the space of states it is clear that the minimal overlap with 
Q can also be reached by an optimal pure state |'0b)('0&|- Now we can choose the following 
bipartite witness 

W b = Q-a b \ifo)(i/> b \, (13) 

and verify that, indeed, for each b , Tr(W&p&) ^ 0 for all p b E Sb and Tr(l / U b |-^ 6 )(^ fe |) = 0, i.e. 
that this operator is a weakly optimal witness for the bipartition {6|fe}. 

Inserting this set of optimal bipartite witnesses into the construction from above using 
[-a h \'fib){fi>b\}+ = 0 yields 

Wqme = Q, (14) 


proving that each weakly optimal multipartite witness can be gained from the construction 
using weakly optimal bipartite witnesses. □ 
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Corollary III.2 Every multipartite entanglement witness Wqme can be constructed using 
the framework summarized by equations flU]) and ([TUP . Furthermore, it is possible to impose 
that, for every bipartition {b\b}, Wb = Aj* <8> 1^[| for some choice of positive map Ab, 
acting on the subsystems in b, and some choice of pure state \iff)- 

Proof: First it is important to notice that, for every state p which is entangled across a 
specific bipartition {6|5}, there exists a positive map Aacting on the subsystems in b, such 
that Ab<g) 1 b [p] 0, i.e. it is detected to be entangled by this positive map (see^—“^). This 
implies in turn that there exists a pure state | ifb) such that (f>b\Ab <8> %[p] It/’fe) < 0, i.e. 

Tr (pA* b ® l b [\ip b )(i} b \]) < 0. (15) 

Through continuity this implies that, for every extremal biseparable state \<f>), there exists 
a weakly optimal witness of the form Wb = A£ ® Ijt[IV' fo)(V'ftl] such that Tr (Wb\4>)((j)\) = 0. 
Now, invoking the Choi-Jamiolkowski map, we can conclude that every possible hyperplane 
intersecting the biseparable set corresponds to a witness derived from a positive map. This 
implies that amongst all the {Wb}b, there is at least one them Wb 0 which satisfies the 
following: for each bipartition {b\b}, Wb = Wb 0 — /3&1, for some fdb , is an optimal bipartite 
witness for p (and /3& 0 = 0). This concludes the proof: the {Wb}b are all obtained from 
shifting W bo = A* o <g) %[|^&o)(^&ol] ■ D 


B. Illustration on one example 


To showcase the strength of the above constructions let us illustrate it with a peculiar 
example. First of all let us make a specific choice for Q that has already proven to work 
well in Ref.-: Given a set of witnesses for every bipartition {Wb}b, we will construct Q by 
finding element wise the matrix of largest common negative matrix entries 


N = ^ | rj) (r/1 min 0, max [iR,e(r]\Wb\ri')] 

r),r]' £{0,l,...,d— 1 } 

and the matrix of smallest common positive matrix entries 


(16) 


P— ^ \r])(r]'\max 0, min [iR.e(p\W b \r]')] 

r/,ri' c{0,l,...,d— 1 } 


( 17 ) 






With these two matrices we can define 


Q:=N + P. 


(18) 


Now for the purpose of elucidating how this construction works in practice let us follow 
it through step by step in an exemplary case. To define our target state a few shorthand 
notations will be useful. Fix a,b,c > 0 and define the following vectors in a three-qutrit 
system 


|Vh) = \A|001) + \/a _1 |110), 

|</> 3 ) = v'a|100) + v'a =T |011), 

|^ 5 ) = Vb\121) + Vb^\ 212 ), 

\ ifo ) = a/S| 220) + \4Fi|002>, 

l^g) = a/c| 022) + Vc^OO), 


\ih) = v^|010) +Va^llOl), 

|Vh> = v / 6|112) + v / F t |221), 

|^ 6 ) = Vb\211) + Vb^\ 122 ), 

|^s) = x/c| 202) + x/cFi 1020), 

|Vho) = 1000) + | 111) + 1222). (19) 


We can then construct a mixed state p as follows 

10 

P = +P(|001)<001| + 1) 


P ■ = 


1=1 
P 


Tr p 


( 20 ) 


By construction this state is PPT across all three cuts. If we choose a = 10 -6 , b = 300 and 
c = 12 x 10~ 3 , then for p > 0.0003 it is also positive under application of the Choi map (or 
its dual) across all three cuts. As it is both PPT and Choi-positive it is fair to say that if 
it is entangled it is only very weakly so. The system size is still small enough to apply the 
map-mixer introduced in the previous section, and to reveal that it is nonetheless indeed 
genuinely multipartite entangled for various values of p > 0.003. 

Hence, this is a good opportunity to demonstrate the power of witness liftings. We can 
use the indecomposable map introduced in Ref.—, choosing c 1 = 1, c 2 = 10 -3 and c 3 = 10 3 , 
revealing entanglement across all three cuts. In fact, even without ever calculating any 
eigenvalues, we can just apply the map on the three subsystems of |Vho) only- That is 


W\ — A?’ C2 ’ C3 <g) I23[|f/ho) (f/ho|]j 
W 2 = Ar 2 ’ C3 ®l 13 [|tM(Vho|], 

W' 3 = A^ ca - cs ® li 2 [|^i 0 )(^ 10 |]. (21) 


9 






Plugging the three resulting witnesses in our construction for Q we get a GME-witness that is 
able to reveal genuine multipartite entanglement in the state for a range of 0 ^ p ^ 0.00069. 
This example illustrates that even bound entangled states which are positive with respect 
to paradigmatic maps can exhibit multipartite entanglement. As the state itself is full rank 
and not symmetric, there is no other known method that could have revealed it to be GME, 
clearly demonstrating the power of positive map mixers and witness liftings. 


C. Finding a multiparticle witness with semidefinite programming 

In Ref.- the multipartite witness Wgme was constructed by starting with bipartite wit¬ 
nesses Wb for each bipartition {£>|6}, and then one possible choice for the operator Q was 
explicitly constructed. While the presented choice works well for many examples, as illus¬ 
trated in the previous section, it is not clear whether it is optimal in the general case. For 
a given state p, however, the optimal multipartite witness Wqme can directly be computed 
as a semi-definite program that is easier to run than the map-mixers. 

For that, consider the following constrained optimization problem 


minimize: 

Ti^pWqme) 

(22) 

subject to: 

V b , Wgme > W b . 

(23) 


This is a semidefinite program, which can be easily and efficiently solved using standard 
numerical techniques. 

There is, in addition, a variation of the presented method for obtaining a multiparticle 
witness from a set of bipartite witnesses. The condition imposed by Eq. (|8]) guarantees 
that, for each b, Tv^oWgme) ^ Tr(cxH4) for any state a, so that Wgme is indeed an 
entanglement witness. However, for being an entanglement witness it suffices that, for each 
6, Tr(abWGME) ^ Tr(oyH4) for any state a b which is separable for the bipartition {b\b}. And 
this is already guaranteed if, for all b , there exists a positive map A b , acting on subsystems 
in b, such that 

A b G> ^ A b (g) % [W&]. (24) 

So, for computing a multipartite witness for a given state one can also consider the semidef- 
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inite program 


minimize: T ic(pWgme) (25) 

subject to: V 6, A b <g) t b [W GME \ ^ A b <8> 1 b [W b \. (26) 

These two semidefinite programs can be used to construct the best witness for a given state 
systematically. They might be useful, if the analytical method from Ref.- does not work. 

It should be noted, however, that the presented formulations are not necessarily the 
best way to construct a multiparticle witness from bipartite witnesses, if one wishes to use 
semidefinite programming. The reason is the following: For each bipartite witness one can 
construct via the Choi isomorphism a positive, but not completely positive map. This map 
detects more states than the original witness. Given these maps one can then evaluate the 
corresponding map relaxations, as described in Section 2. This criterion will be stronger 
than the semidefinite programs presented above. For that reason, we do not discuss detailed 
examples here. 

IV. RELAXATIONS OF SEPARABILITY BEYOND POSITIVE MAPS 

So far, we considered only the approximation of the biseparable sets by super-sets which 
are associated to a positive map. One can, however, also use other separability criteria, 
such as the computable cross norm, aka realignment, (CCNR) criterion^ 1 ^ or symmetric 
extensions^. In the following we explain how the CCNR criterion can be used. 

A. Description of the method 

Let us start by explaining the CCNR criterion. Any bipartite quantum state can be 
expressed via the Schmidt decomposition in operator space as 

Pab = X ® G f . (27) 

k 

where the A*, are positive coefficients, and the G k and Gf are orthogonal observables, i.e. 
they fulfill Tr (GfGf) = Tr (GfGf) = Sij. With this representation one can easily prove 
that the following holds: 

Pab separable =>■ ^ 1. (28) 

k 
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And this necessary condition for separability is known as the CCNR criterion. The criterion 
has the advantage that it detects entanglement in many states where the PPT criterion fails. 
On the other hand, not all two-qubit entangled states can be detected by this test. 

From this structure, one can easily write down entanglement witnesses. Namely, any 
operator of the form 

W = t-'£Gi®G B k (29) 

k 

is an entanglement witness, as it is positive on all states with Y2k Xk ^ 1. This structure of 
the witnesses can be used for constructing witnesses for genuine multiparticle entanglement 
as follows. Consider a witness which has for any bipartition {b\b} the structure 

Wgme = A + G i ( 3 °) 

k 

where G\ and G b k are orthogonal observables on the bipartition {b\b}, and A> ^ 0 is positive 
semidehnite. Clearly, if a state obeys the CCNR criterion for some bipartition, the mean 
value of the witness Wgme will not be negative. Consequently, the witness is also non¬ 
negative on all biseparable states. 


B. Example: The three-qubit GHZ state 

The witnesses from the CCNR criterion are more difficult to handle than the witnesses 
from positive maps. The reason is that no approach via semidehnite programming is possible. 
Moreover, the condition from Eq. (1301) is more difficult to check than the condition in Eq. (JHJ) . 
Nevertheless, we will present an example where known optimal entanglement witnesses have 
this structure. 

Consider first the three-qubit Greenberger-Horne-Zeilinger (GHZ) state | GHZ) = (|000) + 
|111))/x/2- The typical witness for this state is W = 1/2 — \GHZ)(GHZ\. Now, the GHZ 
state can be expressed in terms of its stabilizers as 

\GHZ)(GHZ\ = -(111 + ZZ1 + Z1Z + 1 ZZ + XXX - XYY - YXY - YYX), (31) 
8 

where 1,X, Y, Z represent the Pauli matrices l,a x ,a y ,a z , and tensor product signs have 
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been omitted. After a change of the normalization this can be used to write the witness as 


W = l-2\GHZ) (GHZ | (32) 

[ Z Z1 + 1Z XXX-YY -YYX + XY 1 11 + ZZl 

“ 1 “ |yl Vs + 7! 71 + 7f 71 + 7I Vs J' (33) 

From this representation, it is clear that W is a witness as in Eq. (1301) for the bipartition 
{{1} |{2, 3}}, with Pi = 0. Due to the symmetry, this works for all bipartitions. 

V. ESTIMATING THE PERFORMANCE OF PPT RELAXATIONS IN 
HIGH DIMENSIONS 

In this final section we want to discuss the overall performance of such relaxations in 
multipartite systems, using the paradigmatic partial transpose map. In order to estimate 
the performance of using PPT relaxations to detect randomly chosen multipartite entangled 
states we derive lower-bounds on the fraction of multipartite entangled states, among states 
which are positive under partial transposition across every cut. The latter condition is 
strictly stronger than the relaxation employed in Eq. ([2]), thus providing an upper-bound on 
the fraction of states in 7Z{T b } b (where T b stands for the transposition on the subsystems in 
b) that are also in S ( 2 ). 

Our main result can be summarized as follows: For a fixed number of parties, the ratio 
between the volume of biseparable states and the volume of fully PPT states (as measured by 
either the volume radius or the mean width) scales as 1 /Vd, where d is the local dimension. 

In order to precisely formulate this result, we need to introduce first some of the basic 
notions and definitions that will be employed in the derivations. 

A. Notation and preliminary technical remarks 

We denote by H( C n ) the set of Hermitian operators on C n , on which we define || • ||i r 
as the trace class norm, || • ||hs as the Hilbert-Schmidt norm and || • || op as the operator 
norm. We also denote by || ■ || the Euclidean norm on C n . In what follows, we will be 
mostly interested in the asymptotic regime, when the dimension n tends to infinity. In this 
setting, the letters C, c, c 0 etc. denote numerical constants, independent from any other 
parameters, in particular from n. The value of these constants may change from occurrence 
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to occurrence. When A and B are quantities depending on n, the notation A ~ B means 
that the ratio A/B tends to 1 when n tends to infinity. 

When working with a random variable X , we will use the notation P(„4.(A")) to denote 
the probability of the event A(X), and the notation E (/(AT)) to denote the expectation of 
the function f(X). 

Extra notation, concepts and results from classical convex geometry, which are required 
throughout our proofs, are gathered in Appendix [Aj 

ft may, however, be worth mentioning that whenever we use tools from convex geometry 
in the space "H(C”) (which has real dimension n 2 ) it is tacitly understood that we use 
the Euclidean structure induced by the Hilbert-Schmidt inner product (A, B) = Tt(AB). 
For instance, Definition IA.1I of the volume radius of a convex body K C "H(C n ) becomes, 
denoting by B HS the Hilbert-Schmidt unit ball of "H(C”), 

( vol K \ 1 ^ n 

= (^) - < 34 > 

While Definition IA.21 of its mean width is, denoting by Shs the Hilbert-Schmidt unit sphere 
of "H(C n ) equipped with the uniform probability measure a, 

w(K) — f max Tr(XM) da(X). (35) 

J S H s MgK 

As also mentioned in Appendix [Aj the latter quantity can be re-expressed via Gaussian 
variables, which yields here 

w(K) = — E ( max Tr (GM) ) , (36) 

7n \ M ^ K J 

where G is a matrix from the Gaussian Unitary Ensemble (GUE) on C n and = 
E 11Grills ~„_s. +00 n (see e.g. Ref.—, Chapter 2, for a proof). 

To be fully rigorous, let us make one last comment. All the convex bodies of 7-L(C n ) 
that we shall consider will actually be included in the set D(C n ) of density operators on C n 
(i.e. the set of positive and trace 1 operators on C n ). So we will in fact be working in an 
ambient space of real dimension n 2 — 1 (namely the hyperplane of C n ) composed of trace 
1 elements). This subtlety will not be an issue though, since we will be mostly interested in 
the asymptotic regime n —> +oo. In this setting, the operator that will play for us the role 
of the origin will naturally be the center of mass of D(C n ), i.e. the maximally mixed state 
t/n. 
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Theorem V.l On C n , the volume radius and the mean width of the set V of all quantum 
states satisfy the asymptotic estimates, 

e" 1/4 

vrad('D) ~ —— 

71—>-+00 


n 


w(V) ~ —=. 

rw+oo -v/77. 


(37) 

(38) 

n-r-too yj‘U 

Proof: Eq. (137|) was established in Ref.—. 

Eq. (155|) is a direct consequence of Wigner’s semicircle law (see e.g. Ref.—, Chapter 2, 
for a proof). Indeed, we have by definition 


w(V) = — E I rnaxTr 
7n Vpg© 


G[p- 


i 


n 


) = — E ^maxTrfGpj^j = — E (A m ax(G)) , (39) 

/ ’’in \ PG® J 


where G is a GUE matrix on C n and we denoted by A ma x(G) its largest eigenvalue. The 
claimed result then follows from q n ~ n _j. +0O n and E (A max (G)) ~ n _>. +00 2^/n. □ 

B. Volume estimates 


In the sequel, we shall consider the multipartite system (C d ) 0fe , and slightly adapt and 
generalize the notations introduced in Section [Til We shall denote by S and V the sets of 
states on (C d )® fc which are, respectively, separable and PPT across any bi-partition, and 
by 5(2) and P( 2 ) the sets of states on (C d )® fc which are, respectively, bi-separable and bi- 
PPT. These sets may be more precisely defined in the following way. There are N k = 
2 fc ~ 4 — 1 different bi-partitions of the k subsystems C d . Denoting by {5 1 ,..., S Nk } and by 
{P 1 ,..., V Nk } the sets of states which are, respectively, bi-separable and bi-PPT across one 
of these, we have 

N k N k 

5 = P|5 l and P = p|P\ 


i=\ 


2—1 


'N k 


'N k 


«S( 2 ) = conv [J S 1 


<i=l 


and V( 2 ) = conv I IP’ 

\i=l 


Theorem V.2 There exist positive constants Cd —>d^+oo 1 such that, on (C d )® fc , the volume 
radius and the mean width of the set of states which are PPT across any bi-partition satisfy 

2 k 


w{fP) ^ vrad(P) ^ c d 


#/ 2 ’ 


(40) 


where one may choose c — e 1//4 /4. 
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Proof: The first inequality in Eq. (14(7]) is just by the Urysohn inequality (see Theorem IA.3D . 

To show the second inequality in Eq. (110]) . we will use repeatedly the Milman-Pajor 
inequality (see Theorem IA.6j) and more specifically its Corollary IA.71 We will in fact show 
more precisely that there exist Cd -^-d^+oo 1 such that 

c N k - 1/4 

vrad(P) > c d dk/2 . (41) 

The hrst thing to note is that, denoting by Ti,..., Tjy fc the partial transpositions across the 
A*, different bi-partitions of the k subsystems C d , we have 

v = v n v Tl n • • • n v F n * . (42) 


Now, by Corollary IA.7I applied to the convex body T> C 'H{{C> d )® k ) (which indeed has the 
origin t/d k as center of mass) and to the isometry Tx, we get 


vrad (V D D ri ) ^ 


1 vrad(D) 5 

2 w(V) 


,- 1/4 


,1 .!» C #/2 ’ 


(43) 


the last equivalence being by Theorem IV. 11 We may then conclude recursively that Eq. (144]) 
actually holds. □ 


Theorem V.3 On (C d )® k , the volume radius and the mean width of the set of bi-separable 
states satisfy 

vrad(d>( 2 )) ^ w(S( 2 )) < . ( 44 ) 

where one may choose C = min |6^/ln(l + 2/<5)/(l — 2 5 2 ) 2 : 1/10 < S < 1/4 j and C d ,k = 
y /8 ln(2 )/d k ~\ so that C ^ 11 and Cd,k —^-h-oo 0. 

Proof: The hrst inequality in Eq. (144]) is just by the Urysohn inequality (see Theorem IA.3D . 

To show the second inequality in Eq. (144]) . we start from the following observation: Let S 
be one of the k set of states on (C d )® fc which are separable across a given cut C d : (C d )® fc_1 . 
Then, for each 1 ^ i ^ N^, 

w(S l ) ^ w(S) ^ (45) 

where C = min |6-\/ln(l + 2/<5)/(l — 2h 2 ) 2 : 1/10 < 8 < 1/4}. It relies on the already 
known fact that there exists a universal constant C such that, for any m, n G N with m ^ n, 
the mean width of the set S of separable states on C m ® C n is upper-bounded by C/rriy/n. 
This result was basically proved in Ref.—, Theorem 1, but since specifically stated there in 
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the balanced case m = n only, for vrad(iS) rather than w(S) and without specifying that 
one may choose C = C/ 2, we briefly recall the argument here. 

Let 1/10 < 5 < 1/4 and consider M.$, Af$ 5-nets for || • || within the Euclidean unit 
spheres of C m and C n respectively. Imposing that A 4s, J\fs have minimal cardinality, we 
know by volumetric arguments (see e.g. Ref.—, Lemma 4.16) that \M.s\ ^ (1 + 2/5) 2m and 
\J\fs\ ^ (1 + 2/5) 2 ' 1 . Then, it may be checked that 

conv (S U -S) C conv { ± |s)(a:| ® \y)(y\ : | x) E M s , \y) E Af s j. (46) 

So by Lemma IA.41 we get 


vj(S) V w (conv ( S U —S)) 

1 /2 In (2(1 + 2/5) 2m (l + 2 / 8 ) 2n ) 

^ (1 — 25 2 ) 2 y (mn) 2 

i Vm m + n) In (1 + 2/8) + ln(4) 

(1 — 25 2 ) 2 mn 

3^/111(1 + 2 /S) 

^ (1 — 25 2 ) 2 my / n’ 


(47) 

(48) 

(49) 

(50) 


which is precisely the content of Eq. (J45]). 

Now, we also have that, for each 1 ^ i ^ N^, S * C B tr C Bus- Hence, by Lemma [A.51 
we get 


u;^)) + 2 


( max w(S l ) + 

l^N k 


2 hi (N k )\ 

(d fc ) 2 ) 


C ^81n(2 )k 

d(k+ 1)/2 d k 


c + c d , k 

d( k + 1)/2 ’ 


(51) 


where C d ,k — y/8 ln(2)/d fc_1 . □ 

The conclusion of Theorems IV.21 and IV.31 may be phrased as follows. On a multipartite 
system which is composed of a small number of big subsystems (k fixed and d +oo), 
imposing that a state is PPT across any bi-partition (i.e. the strongest notion of PPT) is 
still, on average, a much less restrictive constraint than imposing that it is bi-separable (i.e. 
the weakest notion of separability). Indeed, the “sizes” of these two sets of states (measured 
by either their volume radii or their mean widths) scale completely differently: the “size” 
of the former is at least of order 1 /d k ^ 2 while the “size” of the latter is at most of order 
1 /G?( fc+1 )/ 2 , hence differing by a factor of order at least \fd. 
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C. A class of fully PPT and GME states 


In Ref.- an explicit class of GME states that were PPT across all cuts was presented. In 
small dimensions it is a hard task to find such examples, but the results from the previous 
section suggest that at least in high dimensions this should be a generic feature of PPT 
states. To emphasize this fact we present an explicit construction of states PPT across all 
cuts and GME with high probability. 

Consider the following random state model on (C rf )® . fix some parameter 0 < a < 1/4 
(independent of d), pick G a traceless GUE matrix on (C d )® fe , and define the “maximally 
mixed + gaussian noise” state on (C d )® fc 

Pa = 5* (* + ff G ) ■ (52) 

Then, typically (i.e. with probability going to 1 as d grows) pc is fully PPT and nevertheless 
GME. 

More quantitatively, we will show that the following result holds. 

Theorem V.4 Let G be a traceless GUE matrix on (C d )® k . Then, the state pc on (C d )® k , 
as defined by Eq. 052D, is fully PPT and not bi-separable with probability greater than 1 — 
exp(— cd k ~ l ), for some universal constant c > 0. 

Theorem IV.dl is a straightforward consequence of Propositions IV751 and [VT6] below. Before 
stating and proving them, let us elude once and for all a slight issue: a GUE matrix on C n 
is the standard Gaussian vector in Tl( C n ), while a traceless GUE matrix on C n is the 
standard Gaussian vector in the hyperplane of TL{ C n ) composed of trace 0 elements. So in 
the asymptotic regime n — y +oo, all the known results on nxn GUE matrices that we shall 
use also hold for traceless nxn GUE matrices (because the ambient spaces of these two 
gaussian vectors have equivalent dimensions in this limit). 


Proposition V.5 Let G be a traceless GUE matrix on (C d ) <glfc . Then, the state pc on 
(C d )® k ; as defined by Eq. (152]) . satisfies 

P (p G ^V)^e~ cdk , (53) 


where c > 0 is a universal constant. 
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Proof: In Ref. 2 -, a deviation inequality is proved for the smallest eigenvalue of a GUE matrix, 
namely: Let G be a GUE matrix on C n and denote by A m i n (G) its smallest eigenvalue. Then, 
for any e > 0, 

p (A min(G) < -(2 + e)y/nj ^ e- ce3/2n , (54) 

where c > 0 is a universal constant. 

Now, observe that G as well as all its partial transpositions G Vi , 1 ^ i ^ N kl are GUE 
matrices on (C d )® k . Hence, Proposition IV.51 follows directly, by choosing for instance e — 1. 
Indeed, by assumption on a, we have 3a < 3/4 < 1, so the probability that pc or any pjd, 
1 ^ i ^ N k , is not positive is less than e~ cd . □ 


Proposition V.6 Let G be a traceless GUE matrix on (C d )® k . Then, the state pc on 
(C d )® k , as defined by Eq. (I52]h satisfies 

P (pg e S (2) ) ^ (55) 

where c > 0 is a universal constant. 


Proof: Our strategy to show Proposition IV.61 is to exhibit a Hermitian M on (c d )® fc which 
is with probability greater than 1 — exp(— cd k ~ l ) a GME witness for the state pc (he. such 
that Tr (p G M) < 0 while Tr (pM) > 0 for any bi-separable state p). 

Note first of all that, on the one hand, 


ETr ( Pg G) 


a 


d 3k / 2 


ETr (G 2 


a 


d-H-oo d? k G 


d 2 k = ad fc/2 . 


(56) 


while on the other hand, by Theorem IV.31 

E /AU Tr(pG) ^ ^vv E[Tr(G2)1 ' /2 ^ JkrP = “ ,t_I|/2 - < 57 ) 

where we used that for G a GUE matrix on C n , E Tr(G 2 ) ~ n _> +00 n 2 and E[Tr(G 2 )] 1 / 2 ~ n ^. +0O 
n (see e.g. Ref.—, Chapter 2, for a proof). 

Let us now show that the functions G i— > Tt(pgG) and G i— y sup pg5(2) Tr(pG) concentrate 
around their respective average values. In that aim, we shall make use of the following 
Gaussian deviation inequality (see e.g. Ref.—, Chapter 2): Assume that / is a function 
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satisfying, for any Gaussian random variables G, H, |/(G) — f(H)\ ^ cf g ,h\\G — H\\hs f° r 
some cjg,h such that E<t g ^ ^ L. Then, for any e > 0, 


P(|/-E./|>e)^e- C0£2 / i2 , (58) 

where Co > 0 is a universal constant. 

Define f:Ge GUE(C n ) ^ Tr(G 2 ), and / s : G G GUE(C n ® C n ) ^ sup peS Tr(pG), for 
any given set of states £ on C n ® C n . We have hrst 

|/(G) - /(if)| = | Tr(GGt) - Tr(HH')\ (59) 

< ||GG t - HH^\\ tr (60) 

^ (||G||//s + ||//||hs)||G — H\\hs, (61) 


where the last inequality is by the triangle inequality, the Cauchy-Schwarz inequality, and 
the invariance of || • ||hs under conjugate transposition, after noticing that GG^ — HH^ = 
G(Gt - i/t) + (G - H)Hl And second 


\h(G)-MH)\ 


sup Tr (pG) — sup Tr (pH) 

pGS pGS 


^sup |Tr(p[G — H])\ 

pGE 


^ IIG — H\\ H s, 


(62) 

(63) 

(64) 

(65) 


where the next to last inequality is because £ is a subset of the trace class unit sphere of 
Hcrmitians on C n <E> C n . 

Hence, the functions / and /s satisfy the hypotheses of the Gaussian deviation inequality 
with L = 2y d k ~ d ^ +00 2 d k and L — 1 respectively. So by the mean estimates (157)1) and (1571) . 
we have that for any 0 < e < 1, 

P (Tr ( Pg G) < (1 - e)ad k/2 ) < exp (-c 0 (ead 2k ) 2 / (2 d k ) 2 ^j = exp ( -c' 0 e 2 d 2k ) , (66) 

P f sup Tr(pG) > (1 + e)Cd (k ~ 1)/2 ) < exp f-c 0 (eGd (fc - 1)/2 ) 2 ) = exp {-^e 2 d k ~ l ) . (67) 

V g5 (2) / V 7 

As a consequence, we have that for any /3 d satisfying 1/2Gd^- 1 )/ 2 ^ /3d ^ 3/2«d fc//2 , the 
Hermitian M = 1 — /S^G on (C d )® fc is a GME witness for p^ with probability greater than 
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1 — exp(— cd k 1 ), where c > 0 is a universal constant. Indeed, choosing e = 1/6 in Eq. (1661) 
and e = 1/2 in Eq. (157)) , we get 

P (tt ( Pg M) > -^] ^ e~ cd2k and P ( sup Tr (pM) < - ) ^ e~ cdk ~\ ( 68 ) 

V 4 / \pes m 4 J 

which concludes the proof. □ 

We may actually say even more on the random state pc defined by Eq. ([52]) . Indeed, 
dehne for all 0 < e < 1 the state p G (e) on ( c d )® k by 

p G (e)=ep G +( l-e)^ = -^(l + e 

Then, what the proof of Proposition IV.61 additionally tells us is that, as long as e ^ 7 /\/d, 
for some constant 7 > 0, p G (e) is with high probability not bi-separable. This means that 
Pg is typically a fully PPT state on (C d )® k which is not bi-separable, and whose random 
robustness of genuinely multipartite entanglement (as defined in Ref.- 2 -) additionally grows 
at least as \fd when d —* + 00 . 

VI. CONCLUSION 

The problem of characterizing genuine multipartite entanglement and biseparability is 
difficult. Therefore, a natural approach lies in the relaxation of the definition of bisepara¬ 
bility: Instead of considering states which are separable with respect to some bipartition, 
one replaces this set by an appropriate superset, e.g. defined by the PPT condition or some 
other positive but not completely positive map. 

In this paper we investigated this approach from several perspectives. First, we estab¬ 
lished how this relaxation approach with positive maps can be evaluated with semidefinite 
programming and how entanglement witnesses can be constructed for this problem. Then, 
we showed that, in principle, also other relaxations, besides those obtained from positive 
maps (e.g. based on the CCNR criterion), are possible. Finally, we studied the accurateness 
of the relaxation approach. We proved rigorous bounds on the volume of the set of bisep- 
arable states as well as on the volume of the set of states which are PPT for any cut. In 
this way, we showed that in the limit of large dimensional multipartite systems, the relax¬ 
ation approach detects only a small fraction of the multiparticle entangled states. It must 
be stressed, however, that this does not mean that the relaxation approach is not fruitful: 


a 

#72 


G 


(69) 
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First, it is a well known fact from the theory of two-particle entanglement that simple en¬ 
tanglement criteria miss most of the states if the dimension of the local spaces increases^ 1 ^. 
Second, from a practical point of view, the relaxation approaches are clearly the best tools 
for characterizing multiparticle entanglement available at the moment^. 

For future research, there are many open questions to address: First, a systematic ap¬ 
proach for the various positive maps besides the transposition would be desirable. Then, 
an approach for characterizing separability classes besides biseparability (e.g. triseparabil¬ 
ity) would be useful. Finally, methods to certify the Schmidt-rank or the dimensionality of 
entanglement^ in high-dimensional systems are needed for current experiments. Investigat¬ 
ing the generic scaling of these quantities could also be of interest. 
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Appendix A: Classical convex geometry 

We work in the Euclidean space R n , where we denote by || • || the Euclidean norm, induced 
by the inner product (•,•). We denote by vol(-) the n-dimensional Lebesgue measure. A 
convex body K C R" is a convex compact set with non-empty interior. 

If u is a vector from the unit Euclidean sphere S n ~ x , the support function of K in the 
direction u is 


h K {u) := max(i, u). (Al) 

x GK 

Note that Iik{u) is the distance from the origin to the hyperplane tangent to K in the 
direction u. 

Two global invariants associated to a convex body K C R n , the volume radius and the 
mean width , play an important role in our proofs. 
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Definition A.l The volume radius of a convex body K C R n is defined as, denoting by B n 
the unit Euclidean ball of R", 

vrad(/ °(^) ■ (A2) 

In words, vrad(A") is the radius of the Euclidean ball with same volume as K. 

Definition A.2 The mean width of a subset K C R n is defined as 

w(K) \= / max(x, u) dcr(w), (A3) 

Js n ~ i X&K 

where dcr(w) is the normalized spherical measure on the unit Euclidean sphere S n ~ l of R n . 
If K is a convex body, we have 

w(K) := f h K (u)dcr(u). (A4) 

J s n ~ l 

The inequality below (see e.g. Ref.—, Corollary 1.4) is a fundamental result which com¬ 
pares the volume radius and the mean width. 

Theorem A.3 (Urysohn inequality) For any convex body K C R n , we have 

vrad(A) ^ w(K). (A5) 


It is convenient to compute the mean width using Gaussian rather than spherical inte¬ 
gration. Let G be a standard Gaussian vector in R 71 , i.e. such that its coordinates, in any 
orthonormal basis, are independent and following a Gaussian distribution with mean 0 and 
variance 1. Setting = E ||G|| ~ y/n, we have, for any compact set K C R n , 

(G,x)^j = 'y n w(K). (A6) 

We also need the two following lemmas, which are incarnations of the familiar “union 
bound”. They appear for example in Ref.—, Chapter 3, under the equivalent formulation 
via suprema of Gaussian processes. 


wa(K) '■= E ( max 

v 2 V x&K 


Lemma A.4 (Bounding the mean width of a polytope) Let Vi,... ,vn be points in 
R n such that u* G A B n for every index 1 ^ i ^ N (where B n denotes the unit Euclidean ball 
ofTU 1 ). Then 


w (conv(vi ,..., v N )) ^ A 


2 In A 


n 


(AT) 
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Lemma A.5 (Bounding the mean width of a union) Let Ad, ... , Abv be convex sets 
in R n such that K\ C A B n for every index 1 ^ i ^ N (where B n denotes the unit Euclidean 
ball of R n ). Then 

( , , , 12 In a\ /a s 

^ 2 I max w(Ki ) + X\ - ) . (A 8 ) 

Finally, we use repeatedly the following result, established in Ref.—, Corollary 3. 

Theorem A.6 (Milman—Pajor inequality) Let K , L be convex bodies in R n with the 
same center of mass. Then 

vrad(A" n L)vrad(A — L) vrad(A")vrad(L), (A9) 

where K — L = {x — y : i 6 A', y e L} stands for the Minkowski sum of the convex bodies 
K and —L. 

Choosing K = —L in Theorem IA. 6 I yields the following corollary. 

Corollary A. 7 If K is a convex body in R n with center of mass at the origin, then 

vrad(AT D - K ) ^ ivrad(A^), (A10) 

and more generally for any orthogonal transformation 9, 

vrad(ATlO(/f))^t ™^) 2 (AH) 

The latter inequality is simply because, on the one hand, vrad( 6 *(A')) = vrad(A), and on 
the other hand, by Theorem IA.31 vrad(AT — 9{K)) ^ w(K — 9(K)) ^ w(K) + w(9(K)) = 
2w(K). 

We typically use Corollary IA.7I in the following way: if AT is a convex body with center 
of mass at the origin which satisfies a “reverse” Urysohn inequality, i.e. vrad(A') ^ aw(K) 
for some constant a, we conclude that the volume radius of AT n 9{K) is comparable to the 
volume radius of K. 
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